Ultracold collisions involving heteronuclear alkali metal dimers 
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We have carried out the first quantum dynamics calculations on ultracold atom-diatom collisions 
in isotopic mixtures. The systems studied are spin-polarized ^Li -I- ''Li^Li, ^Li -I- ®Li2, ®Li -I- 
and ®Li -f ^Li2. Reactive scattering can occur for the first two systems even when the molecules 
are in their ground rovibrational states, but is slower than vibrational relaxation in homonuclear 
systems. Implications for sympathetic cooling of heteronuclear molecules are discussed. 
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There is at present great interest in the properties 
of cold and ultracold heteronuclear molecules. There 
are several reasons for this. First, there are new 
phases of matter predicted for quantum-degenerate dipo- 
lar gases. Secondly, heteronuclear molecules can be 
fermionic whereas neutral homonuclear molecules can- 
not. Dipolar molecules have potential applications in 
quantum computing and in the emerging field of ultra- 
cold chemistry. 

Heteronuclear alkali metal dimers, which can be cre- 
ated from ultracold atomic gas mixtures, are of partic- 
ular interest. Formation of such dimers by photoasso- 
ciation has already been detected in cold atomic traps 
Qi IE IE ^ IE IE Q ' Very recently, magnetically tunable 
Feshbach resonances have been observed experimentally 
in two such systems [E • Efforts are under way to use 
these for production of ultracold molecules. Cold bialkali 
dimers have also been formed on helium droplets (IQj . 

The molecules (heteronuclear or homonuclear) pro- 
duced so far from ultracold atomic gases have been 
formed in high vibrational states. These are prone to un- 
dergo inelastic collisions (quenching) with the remaining 
atoms or other molecules; the resulting recoil produces 
trap loss. In the case of homonuclear molecules, this dif- 
ficulty was overcome by making use of fermionic isotopes 
(^Li and ^°K), and tuning the scattering length to very 
large positive values. Pauli blocking then suppressed the 
inelastic collisions significantly, and long-lived molecular 
Bose-Einstein condensates were produced at the end of 
2003 [HllEIIl- 

Efforts are being made to produce alkali metal dimers 
in low-lying vibrational states, either directly or by trans- 
ferring the population from high-lying states. However, 
we have shown computationally that Pauli blocking does 
not su ppr ess inelastic collisions for low-lying vibrational 
states |l4L ITsf . For such states the quenching rates dom- 
inate the elastic rates below 1 mK even when fermionic 
isotopes are used. To trap homonuclear alkali metal 
dimers with small bond lengths, it will therefore be essen- 



tial to produce them in their ground rovibrational state. 

In the case of heteronuclear molecules the situation 
is even more complicated. Even if heteronuclear dimers 
are produced in their ground rovibrational states, the 
molecules may not be stable against collisions. For ex- 
ample, the spin-polarized reaction 



^Li^Li(w = 0, n 0) + ■^Li ^ ^Li -I- ^Li2(t; = 0, n 



0) 
(1) 

is exothermic by 1.822 K because of the difference be- 
tween the zero-point energies of the two dimers. 

The focus of this Letter is the collision of lithiimr atoms 
with lithium dimers in cold atomic gases containing both 
the ^Li and ^Li isotopes. We have carried out the first full 
quantum dynamics calculations for spin-polarized colli- 
sions involving heteronuclear alkali metal dimers. We 
find that quenching via process , is significantly slower 
than the vibrational relaxation processes studied previ- 
ously in homonuclear collisions '14J|. In the present case 
the elastic cross section dominates the inelastic cross sec- 
tions in a region between 10 /iK and 10 mK. 

Quantum scattering calculations require accurate po- 
tential energy surfaces. We have previously shown that 
the interaction potential for s pin- polarized (quartet) Li 
+ IA2 is highly non-additive [ig, with a well depth 4 
times that of the sum of three Li-Li triplet pair poten- 
tials. In the present work, we use a global quartet Lia 
potential obtained from all-electron coupled-cluster elec- 
tronic structure calculations. This is the same potential 
as was used in our work on homonuclear Li -I- Li2 colli- 
sions and will be described in detail elsewhere 17]. 

The methods we use to carry out quantum dynamics 
calculations on systems of this type have been described 
in our previous work on Na -I- Na2 ^E a-nd else- 
where in the context of thermal reactive scattering 2^ , 
so a brief summary will suffice here. The potential energy 
surface is barrierless, so that it is essential to take reactive 
(atom exchange) collisions into account. The positions of 
the nuclei are described in hyperspherical democratic co- 
ordinates. The configuration space is divided into inner 
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and outer regions, and the boundary between them is 
placed at a distance such that coupUngs due to the resid- 
ual atom-diatom interaction can be neglected outside the 
boundary. In the inner region (hyperradius p < 45 ag in 
the present case) , we obtain the wavefunction for nuclear 
motion by propagating a set of coupled equations in a 
diabatic-by-sector basis that is obtained by diagonalizing 
a fixed-p reference Hamiltonian in a basis set of pseudo- 
hyperspherical harmonics. In the outer region, we use 
the Arthurs-Dalgarno formalism which is based on 
Jacobi coordinates, and compute by inwards integration 
regular and irregular solutions of a radial Schrodinger 
equation which includes the isotropic {R~^) part of the 
interaction. Matching between wavefunctions in the in- 
ner and outer regions yields the scattering S-matrix. 

The present calculations deal with systems containing 
two nuclei that are identical and a third that is different. 
We consider spin-stretched states of all three atoms in- 
volved (states with \Mp\ = F = Fmax), for which the nu- 
clear spin wavefunction is symmetric with respect to any 
exchange of nuclei. However, the electronic wavefunc- 
tion for the quartet state is antisymmetric with respect 
to exchange of nuclei because it depends parametrically 
on nuclear positions. Thus the wavefunction for nuclear 
motion must be antisymmetric with respect to exchange 
of the two identical bosonic nuclei (for fermionic alkali 
metal atoms) and symmetric with respect to exchange 
of the two identical fermionic nuclei (for bosonic alkali 
metal atoms). 

Asymptotically, the hyperspherical functions correlate 
with atom-diatom functions; the diatom functions for 
Hund's case (b), which is appropriate for Li2 ('^S+), 
are labelled with a vibrational quantum number v and 
a mechanical rotational quantum number n; n couples 
with the diatomic electron spin s = 1 to give a resul- 
tant j, but for Li2 the splittings between states of the 
same n but different j are very small and are neglected 
here. In the "^E^ state, only even n is allowed for ''Li2 
and only odd n for ^Li2. In the outer region, we in- 
clude rovibrational states with v — 0,1,..., 7 with all 
rotational levels up to Umax = 32, 29, 26, 23, 20, 17, 13, 7 
for ^Li^Li, with all even rotational levels up to n^ax = 
32, 30, 28, 24, 22, 18, 14, 10 for '^Li2, and with all odd ro- 
tational levels up to Umax = 31, 27, 25, 21, 19, 15, 11, 7 for 
^Li2. Because of the reduced symmetry, the present cal- 
culations are substantially more expensive than our pre- 
vious ones on collisions involving 3 identical atoms, and 
we have restricted them to total angular momentum and 
parity — and 1". The basis sets of pseudohy- 
perspherical harmonics involved between 3660 and 6488 
functions. 

There are four collision systems of interest to ^Li/^Li 
mixtures: 



•^Li 
^Li 



^Li2 
'Li2. 



(4) 
(5) 



We have calculated the low-energy elastic, inelastic and 
reactive cross sections <7{E) for all four of these colli- 
sion systems, considering all energetically allowed pro- 
cesses for a variety of initial diatomic rovibrational states 
(f , n). The present paper will focus on collisions with the 
diatomic molecule initially in its lowest allowed rovibra- 
tional state. Details of the results for other initial states 
may be found in ref. i22i . 

At very low energies (up to about 1 mK), the cross sec- 
tions are dominated by collisions with orbital angular mo- 
mentum I — 0, for which there are no centrifugal barriers. 
The s-wave cross sections for ^Li + ^Li^Li(w = 0,rt = 0) 
are shown as a function of collision energy E in Fig. 
n There is an energetically allowed reactive pathway to 
form ^Li2(u = 0, n = 0) + ^Li; at collision energies above 
0.477 K the v = 0,n = 2 state of ^Li2 is also accessible. 




FIG. 1: Elastic and reactive s-wave cross sections for 
''Li^Li(u = 0,n = 0) 



Li + 



^Li + ^Li^Li 
^Li + ^Li^Li 



(2) 
(3) 



Scattering calculations for Li + Li2 are converged up 
to about 10~^ K using just a single partial wave. At 
higher energies, other partial waves start to contribute 
significantly. It may be seen in Fig. ^ that the reac- 
tive cross section for ''Li -I- ^Li''Li(w = 0, n = 0) at 
10~'* K is 6.7 X 10"^'' cm^, corresponding to a limit- 
ing low-temperature rate coefficient fcioss = 4.1 x 10^^^ 
cm'' s^^. This is about 2 orders of magnitude smaller 
than obtained previously for ^Li -I- ^Li2(z; > 0) and 
■"Li + ''Li2(w > 0) Q; the difference can be attributed 
to the smaller kinetic energy release and single product 
rovibrational channel in the present case (and is present 
only for the initial state v ~ 0,n — 0). As a result, 
elastic scattering dominates reactive scattering for ^Li -|- 
®Li^Li(w = 0,n = 0) at collision energies between 10 /iK 
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and 10 mK. (The upper limit is uncertain because the 
partial wave sum is not converged at higher energies). 

The existence of a low-energy closed product channel 
increases the density of low-energy Feshbach resonances. 
These are most conveniently detected by plotting the S- 
matrix eigenphase sum <;{E) |23l| as a function of energy: 
at a resonance, <;{E) increases sharply by tt, superim- 
posed on a falling background. The eigenphase sum for 
^Li + ^Li^Li is shown as a function of energy above the 
?; = 0, n = threshold in Fig. |21 There are clear isolated 
resonances at 225 and 470 mK and a pair of overlap- 
ping resonances at about 630 and 725 mK. The cusp at 
825 mK corresponds to the opening of the {v = 0,j = 1) 
reactant channel. Small changes in the potential energy 
surface might result in significant changes in the positions 
of individual resonances, but the density of resonances 
obtained here should be reahstic. However, it should be 
remembered that the present calculations neglect nuclear 
spin and magnetic fields, which produce additional reso- 
nance structures and ways to tune them. 

2 I . . . , ^ . . 1 . . . , ^ . . 1 . ^ 1 
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FIG. 2: Eigenphase sum for ^Li + ^Li^Li scattering (j" — 
0+). 



For ^Li -I- ^Li^Li(w — 0,n — 0), the reactive channel 
to form ^Li2 is closed at energies below 0.825 K and only 
elastic scattering is possible. The low-energy limit of the 
cross section is 4.7 x 10^^^ cm^. 

For ^Li2, the lowest allowed rotational level is n = 1. 
For ''Li -|- ^Li2(u = 0,n = 1), the lowest partial wave to 
include an Z = incoming channel is therefore = 1". 
The corresponding cross sections are shown as a function 
of collision energy E in Fig. O The species can either col- 
lide elastically or react exothermically to form ^Li^Li + 
^Li with the diatom in its (0,0), (0,1) and (0,2) states. In 
this case the reactive cross section at 10~* K is 7.1 x 10"^'^ 
cm^, which corresponds to a limiting low-temperature 



rate coefficient fcioss = 4.4 x 10""'^^ cm^ s~^. This is con- 
siderably larger than for ^Li -I- ^Li''Li('i; = 0,n = 0), 
though the partial cross section to form the {v ^ Q,n ~ 0) 
product is comparable; the larger value for the total reac- 
tive cross section reflects the larger kinetic energy release 
and increased number of open product channels for ^Li 
+ 6Li2(w = 0,n= 1). 

For ^Li + ^Li2(f = 0,n = 0) the reactive channel to 
form ^Li^Li is closed at energies below 1.822 K and only 
elastic scattering is possible. The low-energy limit of the 
cross section is 1.3 x 10~^^ cm^. 
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FIG. 3: Elastic and reactive cross sections for ^Li -|- ^'Li2{v = 
0,n = 1) with J" = 1". 

The results obtained here have important implications 
for experiments aimed at producing a quantum gas of 
^Li^Li in its ground rovibronic state in an ultracold mix- 
ture of ^Li and ^Li atoms. In order to stabilize the 
molecular cloud against two-body trap losses induced by 
the reactive process the remaining atomic "^Li must 
be removed quickly after ground-state molecule produc- 
tion, so that just the two-species fermionic mixture of 
^Li^Li(w = 0, n = 0) molecules and ^Li atoms is left in 
the trap. One could separate the ^Li atomic cloud from 
the molecular cloud as well, but it might be advantageous 
to keep them both in the trap. It has been observed that 
in a one-component dilute Fermi gas the s-wave colli- 
sions are strongly suppressed because of Pauli blocking 
[i^ l . This will make evaporative cooling of the molecular 
cloud very inefficient. However, to achieve further cool- 
ing of the ^Li^Li{vi = 0, = 0) molecules, a sympathetic 
cooling scheme might be developed. Sympathetic cool- 
ing of fermionic molecules can be approached from two 
directions: interstate or interspecies sympathetic cooling. 

Interstate sympathetic cooling is based on s-wave col- 
lisions between particles of the same species but in two 
different states. It was previously used, for example, to 
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produce two overlapping Bose-Einstein rubidium conden- 
sates (25| and to create a Fermi sea from potassium atoms 
. Unfortunately, for ^Li'^Li there are various inelastic 
processes resulting from reactive molecule + molecule s- 
wave collisions that are likely to produce significant trap 
losses. It will therefore be better to keep the molecules 
in the same state, where Pauli blocking suppresses such 
processes. 

Interspecies sympathetic cooling makes use of s-wave 
collisions between two different species. It has been used 
to produce a Bose-Einstein condensate of atomic potas- 
sium and to create a two-species mixture of quantum- 
degenerate Bose and Fermi gases . In the present case 
we can distinguish two interspecies sympathetic cooling 
schemes. 

The first scheme uses the fact that ultracold collisions 
between spin-stretched ^Li''Li(w = 0, n = 0) and ^Li can 
result only in elastic scattering, because the correspond- 
ing inelastic and reactive processes are energetically for- 
bidden. Thus no undesirable quenching-induced trap loss 
occurs. The elastic collisions will then lead to sympa- 
thetic cooling between the atomic and molecular clouds. 
This could be done in the trap where the molecules were 
produced, or alternatively the ^Li^Li molecules could be 
transported into a reservoir of colder ^Li atoms. This 
scheme is similar to interstate sympathetic cooling for 
fermions in the sense that neither of the components can 
be cooled evaporatively and the cooling is done only by 
interspecies collisions. 

If the efficiency of the first scheme is not high enough 
to achieve quantum degeneracy in the molecular cloud, 
then an alternative will need to be adopted. One possibil- 
ity is to transport the molecular cloud into a reservoir of 
bosonic atoms such as ^^Rb. Spin-polarized ultracold col- 
lisions between ^Li^Li(ti = 0, n = 0) and an alkali metal 
atom heavier than Li can lead only to elastic scattering, 
because the reactive quenching channels are not energet- 
ically accessible. The bosonic component can then be 
cooled evaporatively. Note however that this will work 
only if both species are in spin-stretched states, as col- 
lisions involving other states could produce reaction to 
form singlet states of molecules such as LiRb. 
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